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Electrostatic mapping of nuclear pairing 
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The traditional nuclear pairing problem is shown to be in one-to-one correspondence with a clas¬ 
sical electrostatic problem in two dimensions. We make use of this analogy in a series of calculations 
in the Tin region, showing that the extremely rich phenomenology that appears in this classical 
problem can provide interesting new insights into nuclear superconductivity. 
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Pairing is a pervasive feature in nuclear structure. Per¬ 
haps its most dramatic manifestation is nuclear super¬ 
conductivity, the analogue of the more familiar super¬ 
conductivity that arises due to the pairing of electrons in 
a solid. In a nucleus, however, unlike in a solid, only a 
small number of active particles can correlate under the 
influence of a pairing force. As a consequence, all fea¬ 
tures of the transition to nuclear superconductivity are 
softened and it is extremely difficult to see manifestations 
of the superconducting phase transition. In this work, 
we develop an exact mapping between the nuclear pair¬ 
ing problem and a classical electrostatic problem in two 
dimensions. This mapping permits us to recast the prob¬ 
lem of nuclear pairing in a completely new setting, one 
in which the phase transition to superconductivity can 
be vividly seen. As we will see, the rich phenomenology 
that appears in the classical problem under appropriate 
conditions can provide important new insight into the 
phenomenon of superconductivity in nuclei. 

The key to this new view of nuclear pairing comes from 
the fact that the traditional pairing problem is exactly 
solvable. Its solution was first presented by Richardson 
[0 |J in the 1960s in a series of papers. Though recent 
work has shown how to generalize Richardson’s solution 
to other systems governed by pairing §, and there have 
been recent applications of these generalized models to 
boson systems 0 |t) , we will restrict ourselves here to the 
original fermion problem, for which the hamiltonian takes 
the familiar form 
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Note that the pairing hamiltonian (ill) includes a single¬ 
particle energy term that splits the set of active levels. 

As discussed in ref. ||, the eigenstates and eigenvalues 
of the pairing problem can be obtained by solving a set 
of simultaneous eigenvalue equations 


where the Rj are a set of mutually commuting operators 
that depend on a set of parameters {?y j,g}. The pairing 
hamiltonian ( 0 ) corresponds to the linear combination 
H = 2 J2j e jRj + cfe, obtained by choosing the rjj to be 
the single-particle energies. It is worthwhile to note here 
that other exactly-solvable pairing hamiltonians can be 
obtained by making use of the freedom in the parameters 
T] 0,01. The form of the R operators was first found in 
Ref. [g[ , where integrability of the pairing hamiltonian ( 0 ) 
was demonstrated. The eigenvectors of the R operators 
are equivalent to those originally proposed by Richardson 
to diagonalize the hamiltonian ( 0 ) , viz: 
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where M is the number of collective pairs and v = y~b v.i 
is the number of unpaired particles, so that the total 
number of particles is N = 2 M + v. [The ground state 
corresponds to v = 0 (y = 1) for even (odd) systems]. 
The M unknowns e a , which are referred to as the pair 
energies , are determined by solving the set of equations 
(!)• 

Solving ( 0 ) leads to two sets of equations. One is a set 
of equations for the unknown pair energies e a (called the 
Richardson equations), 
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where kj = ^ ^ and O, = j + 1/2. The quantity — 2kj 

plays the role of an effective pair degeneracy giving the 
maximum number of pairs, allowed by the Pauli princi¬ 
ple, that the single-particle orbit j can acommmodate. 
The second gives the eigenvalues of the R operators, 
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The pairing eigenvalues are given by E = e °n making 
clear why the e a are called the pair energies. 
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We now discuss how to establish an exact electrostatic 
analogy for the pairing problem, building on ideas first 
published by Richardson |J. To do so, we define an en¬ 
ergy functional U 

U = — i e a + 2 kj£j I — kj In \2ej — e a \ 
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A similar energy functional has been recently derived 
from Conformal Field Theory It is straightforward 
to verify that the Richardson equation (j|) can be ob¬ 
tained by taking derivatives of U with respect to the pair 
energies e a . Likewise the eigenvalues of the R operators 
(jh) in units of g can be obtained from V by taking deriva¬ 
tives with respect to twice the single-particle energies 2e*. 

In searching for a physical meaning to the energy U , 
we remind the reader that the Coulomb potential in 2D 
due to the presence of a unit charge at the origin is given 
by the solution of the Poisson equation 

Av (r) = — 27 t< 5 (r) . (8) 

,/From this, it is easy to confirm that the Coulomb inter¬ 
action between two point particles is 

v (ri,r 2 ) = -gig 2 ln|ri - r 2 | , (9) 

where qi is the charge and r, the position of particle i. 

Returning to the energy functional (^J), we now recog¬ 
nize it as describing a 2D classical electrostatic system 
of L fixed charges (L is the number of orbitals in the va¬ 
lence space) at positions 2e, and with charges ki (these 
charges are in general negative for small values of the 
orbital seniority vf), and M free charges located at po¬ 
sitions e a and with positive unit charge. The real axis 
is mapped onto the vertical axis and the imaginary axis 
onto the horizontal axis. Besides the Coulomb interac¬ 
tion between all charges (the second, third and fourth 
terms in eq. (|t|)), there is a uniform electric field in the 
vertical direction (the first term) whose intensity is in¬ 
versely proportional to twice the pairing strength, 2 g. 
Since the fixed charges are related to the single-particle 
orbitals, we will call them orbitons. The free charges 
associated with the pair energies will be called pairons. 


Table I: Correspondence between properties of the 
quantum pairing model and those of the classical 
two-dimensional electrostatic problem. 


Pairing Model 

Electrostatic 

Single-particle energy: £i 

Orbiton position: 2£j 

Orbital pair degeneracy: f li 

Orbiton charge: ki 

Pair energy: e a 

Pairon position: e a 

Pairing strength: g 

Electric field intensity: 


Table I summarizes the translation of the quantum 
pairing model to the classical electrostatic problem. 

Solving the Richardson equations for the pair energies 
e a is thus completely equivalent to finding the stationary 
solutions for the pairon positions in the classical electro¬ 
static problem. 

Since the orbiton positions are given by the real single¬ 
particle energies, they are fixed to lie on the vertical axis. 
The pairon positions are not of necessity constrained to 
the vertical axis, but rather must be reflection symmet¬ 
ric around the vertical axis. This reflection symmetry 
property can be readily seen by performing complex con¬ 
jugation on the Coulomb energy functional (!7f). As a 
consequence, a pairon must either lie on the the vertical 
axis (real pair energies) or must be part of a mirror pair 
(complex pair energies). The various stationary pairon 
configurations can be readily traced back to the weakly 
interacting system (g —> 0). In this limit, the pairons are 
distributed around the orbitons forming artificial atoms. 
As mentioned above, the number of pairons surround¬ 
ing the ith orbiton cannot exceed — 2 ki, which for the 
ground state (where v = 0) is fThus, for small g , 
the ground-state configuration corresponds to distribut¬ 
ing the pairons around the lowest-position orbitons con¬ 
sistent with this Pauli constraint. We then let the system 
evolve gradually with increasing g until we reach its phys¬ 
ical value. The excited states can be constructed start¬ 
ing from a different initial pairon configuration and/or 
by breaking pairs (increasing the seniority v). 

As an example, we now show some results for the 
ground states of the two semi-magic nuclei 114 Sn and 
the 116 Sn. Details on the solution of the Richardson 
equations (||) can be found in ref. ||. Table II shows the 
orbiton positions and charges, which are the same for 
both nuclei under the usual assumption that the single¬ 
neutron energies do not change with neutron number. 

Table II: Positions and charges of the orbitons 
appropriate for a description of the Sn isotopes. 


Orbiton 

Position 

Charge 

^5/2 

0.0 

-1.5 

97/2 

0.44 

- 2.0 

s l/2 

3.80 

-0.5 

^3/2 

4.40 

- 1.0 

^11/2 

5.60 

-3.0 


The nucleus 114 Sn has 14 valence neutrons and thus in 
our classical electrostatic analogy seven pairons. In the 
weak coupling limit, the lowest configuration has three 
pairons close to the d 5 /2 orbiton and the other four close 
to the < 77/2 orbiton. Because of the reflection symmetry 
property, the four pairons close to the < 77/2 orbiton form 
two mirror pairs, whereas the three that are close to the 
d 5 / 2 orbiton form one mirror pair and an odd pairon con¬ 
strained to the vertical axis. 
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FIG. 1. Two-dimensional representation of the pairon posi¬ 
tions in the 114 SVi for three selected values of g. The orbitons 
are represented by open circles; only the lowest two, the d 5 / 2 
and P 7 / 2 , are shown at the positions dictated by Table II. 

In figure 1 we show the pairon positions in the 2D plane 
for three selected values of g. [The physical value of g is 
approximately —21.7 MeV/A = —0.190 MeV.} As seen 
in the figure, for small values of g the pairons are in¬ 
deed organized around the two lowest orbitons and form 
two atoms. The other three orbitons correspond to levels 
that are fairly high in energy and thus lie outside the fig¬ 
ure. The fact that the pairons tend to stay in positions 
below the corresponding orbitons is due to the interplay 
between the orbiton attraction and the external electric 
field which points downwards for attractive pairing. As 
g increases, the atoms expand to reduce the Coulomb in¬ 
teraction energy thereby balancing the reduction in the 
electric field and there appears a transition from the two 
isolated atoms to a cluster at around g ~ 0.08. [Note: 
In the figure, each pairon is connected by a line to the 
one that is nearest to it.] We claim that this geometrical 
transition from atoms to clusters in the classical problem 
is a direct reflection of the superconducting transition in 
the quantum problem. The transition from normal to 
superconducting, which is so difficult to see in the exact 
quantum results because of the finiteness of the Hilbert 
space, shows up very clearly and in a highly pictorial 
fashion through the electrostatic analogy. 

In Figure 2, we show the corresponding results for the 
nucleus 116 Sn. At fairly weak coupling (g = 0.06 in Fig. 
2A), seven of the pairons are organized in two atoms, as 
in 114 SYi, while the eighth lies close to the next orbiton, 
the S\/ 2 - This last pairon is then constrained to move 
downwards along the vertical axis as g increases. 
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- 2-10 1 2 



FIG. 2. Two-dimensional representation of the pairon posi¬ 
tions in 116 Sn for three selected values of g before the collapse 
(Panel A) and three values of g after the collapse (Panel B). 
The orbitons are represented by open circles at the positions 
given in Table II. In Panel A, the energy scale extends high 
enough to include four orbitons, the d 5 / 2 , <?7/2, S1/2 and ^3/2- 
In Panel B, only the d 5 / 2 and <? 7 / 2 are within scale. 

The lowest two atoms display a transition to a clus¬ 
ter, already clearly evident by g = 0.12. At this point, 
however, the last pairon is not yet correlated with the 
cluster, but rather still approaching it from above. Since 
it is constrained to be on the vertical axis, there is a crit¬ 
ical value of g at which it crosses the g 7 / 2 orbiton, which 
for this problem occurs for g c = 0.19044, very near the 
physical value of 0.187. At this point, a dramatic phe¬ 
nomenon takes place, as illustrated in Figure 3. To cancel 
the logarithmic divergence in the energy due to the on¬ 
site interaction between the pairon and the < 77/ 2 orbiton, 
all other pairons collapse onto the two lowest orbitons. 
After the collapse, a new expansion takes place with in¬ 
creasing g (Fig. 2B). Now, however, the two odd pairons 
associated with the d 5 / 2 and the Si / 2 orbitons form a mir¬ 
ror pair. All eight pairons, now in four mirror pairs, then 
expand in the 2D plane forming a cluster around the two 
lowest orbitons. At this point and beyond, all pairs are 
collectively involved in the quantum superconductivity. 

The electrostatic analogy also enables us to understand 
the interesting pattern of occupation probabilities for 
116 Sn illustrated in the lower panel of Fig. 4. They were 
calculated for the exact solution using the Hellmann- 
Feynman theorem and taking derivatives of the pairing 
hamiltonian with respect to the single-particle ener- 
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HF+BCS codes, and in this way improve their accuracy. 
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FIG. 3. Real (Upper) and imaginary (Lower) parts of the 
pair energies for 116 Sn as a function of g. Results are shown 
only for those with Im(e a ) > 0 . 


FIG. 4. Calculated occupation probabilities for 114 Sn (Up¬ 
per) and 116 5’n (Lower) as a function of g. 


can be found in |j|Q]). Note that as g increases, the occu¬ 
pation of the s 1/2 orbit behaves in a very dramatic fash¬ 
ion, first decreasing precipitously, then stabilizing and 
rising together with those of the uppermost d 3 /2 and 
/in /2 orbits. The precipitous drop is associated with 
the rapid movement of the eighth pairon along the ver¬ 
tical axis, from the S \/2 orbiton, with which it was ini¬ 
tially connected, towards the 37/2 orbiton. As this is tak¬ 
ing place, the other seven orbitons behave exactly as for 
114 SVi, as do the corresponding occupation probabilities 
(see the upper panel of Fig. 4). Once the collapse takes 
place, the eighth pairon suddenly becomes part of the 
cluster. Correspondingly, from this point on all 16 nucle¬ 
ons participate in the superconductivity. For larger val¬ 
ues of g , the S 1/2 orbit behaves much like the two higher 
ones, gradually increasing its occupation as the collectiv¬ 
ity becomes further enhanced. In the large-g limit, all 
occupation probabilities converge to a value of 1 / 2 . 

For moderate dimensions, like the examples discussed 
in this letter, the pairing hamiltonian can be exactly di¬ 
agonalized in a quasispin basis. The Richardson solution 
discussed here has two particular features not present in 
other approaches. On the one hand, since the number 
of Richardson equations (and equivalently the number of 
unknowns) is the number of fermion pairs, the method 
can be readily implemented for extremely large model 
spaces involving many major shells. Thus, for example, 
the exact solution could replace the BCS part of Skyrme 


Equally important is that the Richardson solution leads 
naturally to an electrostatic analogy, and this provides 
useful pictorial insight into pairing phenomena. As we 
have seen in this work, it makes clear that the supercon¬ 
ducting state is realized when collective (Cooper) pairs 
are developed which involve the cooperative participa¬ 
tion of all active orbits with all connection to individ¬ 
ual orbits lost. This phenomenon is exhibited in the 
electrostatic picture as the formation of a pairon clus¬ 
ter. The extensions proposed above in the context of 
Skyrme HF+BCS calculations would permit the use of 
the electrostatic analogy to obtain insight into the pair¬ 
ing properties of more complex nuclear systems. 

In summary, we have developed in this work an exact 
mapping of the nuclear pairing model onto the classi¬ 
cal two-dimensional electrostatic Coulomb problem. The 
classical systems that arise by analogy with the Tin re¬ 
gion were studied, revealing a rich phenomenology with 
important implications for our understanding of nuclear 
superconductivity. Finally, though we focused here on 
applications to nuclei, the general ideas and methods we 
presented can be applied to any fermion system governed 
by pairing correlations. 

Acknowledgments This work was supported in part 
by the Spanish DGES under grant # BFM2000-1320- 
C02-02, by NATO under grant # PST.CLG.977000 and 
by the National Science Foundation under grant # PHY- 
9970749. 


4 
















[1] R.W. Richardson, Phys. Lett. 3, 277 (1963). 

[2] R.W. Richardson and N. Sherman, Nucl. Phys. 52, 221 
(1964). 

[3] R.W. Richardson, J. Math. Phys. 6, 1034 (1965). 

[4] R.W. Richardson, Phys. Rev. 141, 949 (1966). 


[51 J. Dukelsky, C. Esebbag and P. Schuck, Phys. Rev. Lett. 
87, 066403 (2001). 

[6] J. Dukelsky and P. Schuck, Phys. Rev. Lett. 86, 4207 

( 2001 ). 

[7] J. Dukelsky and S. Pittel, Phys. Rev. Lett. 86, 4791 

( 2001 ). 

[8] M.C. Cambiaggio, A.M.F. Rivas and M. Saraceno, Nucl. 
Phys. A 424, 157 (1997). 

[9] R.W. Richardson, J. Math. Phys. 18, 1802 (1977). 

[10] G. Sierra, Nucl. Phys. B572, 517 (2000). 


5 



